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Abstract. We investigate creation, manipulation, and steering of entanglement 
in spin chains from the viewpoint of quantum communication between distant 
parties. We demonstrate how global parametric driving of the spin-spin coupling 
and/or local time-dependent Zeeman fields produce a large amount of entangle- 
ment between the first and the last spin of the chain. This occurs whenever the 
driving frequency meets a resonance condition, identified as "entanglement reso- 
nance" . Our approach marks a promising step towards an efficient quantum state 
transfer or teleportation in solid state system. Following the reasoning of Zueco 
et al. [1], we propose generation and routing of multipartite entangled states by 
use of symmetric tree- like structures of spin chains. Furthermore, we study the ef- 
fect of decoherence on the resulting spin entanglement between the corresponding 
terminal spins. 

1 Introduction 

Interacting spins are ubiquitous in physics, whether it be in magnetic molecules [2, cold 
atoms [3], Penning traps 0] or Josephson-j unction arrays [S]. Since years, the equilibrium and 
non-equilibrium properties of interacting spins have been extensively studied. Phenomena such 
as quantum and classical phase transitions and dynamical localization [6] are well-known in 
the literature. The Heisenberg, Ising, and XY coupling models became paradigmatic models 
in condensed matter physics and statistical physics [7]. 

Needless to say, the dynamics of interacting spins is of quantum nature. The interest in 
those systems has been reinforced along with new prospects in quantum technologies, which 
are expected to exceed dramatically the capabilities of their classical analogs. Typical examples 
are quantum metrology [8] or quantum information processing [9|, with applications such as 
simulators or universal quantum computers. With the perspective of quantum computing in 
mind, spin chains are natural connectors of information between the different parts of the 
"quantum hardware" and the readout device p^HT2] . 

A key ingredient for all these applications is the entanglement between different spins in 
an array. Entanglement has no classical analogue and marks the gap between classical and 
quantum technology. Coming back to quantum communication protocols, the entanglement 
shared by the emitter and receptor of the information ensures a communication superior to 
what would be possible with classical protocols [12] . Consequently, the ability to fully control 



^ e-mail: hanggiOphysik.uni-augsburg.de 



2 



Will be inserted by the editor 



entanglement creation and manipulation is a fundamental requirement. One way to achieve this 
is by using time dependent fields. The benefits of external ac-driving has been recently reported 
in Refs. [T1 [T5HT5] . However, from a technological point of view, it appears rather difficult to 
access each spin separately by applying local magnetic fields. 

In this work, we review and extend our recent results [T1[T3] for entanglement control via 
global external fields acting on all spins in an array. Concretely, we discuss entanglement creation 
between both ends of the chain, which turns out to be optimal at particular driving resonances. 
Once entanglement is created locally, neighboring spins interact such that the spins at the 
ends of the chain efficiently communicate with each other and entanglement routing through 
the chain remains the main task. Here, we extend former studies on entanglement creation via 
modulation of the spin-spin coupling |13j to the case where modulation of the local splittings 
is used instead. 

In the second part of this work, we report our idea of quantum information routing by means 
of ac-fields. In detail, we propose a way to produce multipartite entanglement between many 
distant parties just by joining several "quantum routers" in "beam-splitter" mode. Finally, we 
provide studies on the validity of both driving protocols discussed above under dccohcrcnce. 

The paper is organized as follows. In Sec. [2] we introduce the model and quantify the 
entanglement by means of the concurrence between the outermost spins of the chain. Section [3] 
is devoted to entanglement resonance, while in Sec. S] we introduce the idea of routing and 
entanglement distribution in quantum networks. In Sec. [5] we discuss a possible implementation 
in Penning traps. 



2 Driven spin chains and entanglement computation 

In this work we describe a chain of N spins by the anisotropic XY-mode\ with time-dependent 
parameters using the Hamiltonian |I9j {H = 1) 

^ = 2 E ^"(^x + ^ E [(1 + + (1 - ■ (1) 

n—l n—1 

Both the local energy splittings hi{t) and the interaction strengths J{t) are assumed to depend 
through modulated global fields acting on all spins at the same time. Nevertheless we allow for 
inhomogencities of the fields. A non-zero value of the parameter 7 refers to an homogeneous 
anisotropy. For possible applications in quantum communication the spin chain must have 
open ends. Both the isotropic limit of the model 7 = and the thermodynamical limit N — )■ 
00 can be solved analytically in the time-independent case by means of the Jordan- Wigner 
transformation [TnidD]. On the contrary we are interested in both finite-size effects and time- 
dependent dynamics. We consider the fields to be composed of dc contribution and one harmonic 
component, such that 

hn(t) = e„(/io + hi sinojt) , J{t) = Jq + Ji sin tot. (2) 

For later convenience we move to an interaction picture defined hy X = UQ{t)XUo{t), where 
Uo{t) ~ cxp[— i (y3„(t)crfj and (fn{t) ~ /i„(<)/2. Introducing the ladder operators = 
i((T^ ± ia^), the Hamiltonian ([1]) gets replaced by the interaction-picture Hamiltonian 

H = -{Jo + Ji sin(^t)) ^ ka-^,e'^"(*) + ^-a+^^e-^^-W 

^ n=l ^ (3) 

-f 7(e'^"(*V+a++i +e-^"(*)a-a-+i)" , 

where Z\„ = <^n — ^n+i and i7„ = ipn + fn+i - The first two terms in the sum of Hamiltonian ([5]) 
swap the states of spins n and n-l-I, whereas the latter two terms create and destroy excitations, 
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respectively. As it will turn out, the swapping terms are responsible for entanglement transfer 
while the latter pair of terms constitutes the origin of entanglement generation. 

In the following, we are interested in the bipartite entanglement between two spins in the 
chain. For this purpose one needs to obtain the reduced density matrix of spins j and k. It 
results from the full density matrix of the chain gtot by tracing out all other degrees of freedom, 

Qjkit) = Tr„^j-,fc [£»tot(i)] , j,k , (4) 

where Tr„^j.fe denotes the partial trace with respect to all spins n = 1, 2, ... TV but spins j and 
k. The entanglement shared by spins j and k is then given by the concurrence |21| 

Cjk = max{Ai - A2 - A3 - A4, 0} . (5) 

The X[y are the ordered square roots of the eigenvalues of the matrix Qjk{o'^crk)Qjki'^^'^k)- 

3 Resonances and entanglement creation 

The time dependence of the coefficients in the Hamiltonian ([T]) generally leads to a complex 
dynamics, in addition to the nontrivial dynamics for the time-independent case. Therefore, we 
confine ourselves to scenarios in which both amplitudes Ji and hi from Eq. ^ are smaller 
than the dc field strength Hq. We focus on the spectral response of the chain to periodically 
modulated coupling. All other parameters have arbitrary but fixed values. For simplicity, we 
assume in this section a homogeneous chain, e„ = 1 for all n. 

In the protocol, the spins are initially uncoupled, i.e. J{t) = for t < 0, and cooled down 
to the fully-aligned separable state 

1^(^=0)) = 10000...), (6) 

which is the ground state of the Hamiltonian ([ij with J = 0. At t = 0, we switch on both 
the coupling J(t) and the sinusoidal driving h{t). This procedure is termed a quantum quench 
and constitutes a typical protocol for entanglement generation [35] . Even though the interaction 
strength is usually considered as constant, we take into account the possibility of time-dependent 
coefficients Ji{t) > for t > 0. In this sense we are introducing time-periodic quantum quench 
dynamics. 

In the following, we raise the question of how much entanglement between the ends, mea- 
sured in terms of the concurrence, may be created of the chain. Such large-distance entan- 
glement could be used for typical quantum communication protocols, such as teleportation or 
state transfer [12] . 

3.1 Time-dependent interaction 

At first, we only consider the interaction to depend on time, i.e. hi = and Ji 7^ 0. We 
investigate the time evolution of the concurrence between spin zero {n = 0) and the last 
spin (n = TV) for different driving frequencies w^, and the parameters 7, Ji and Jq by direct 
numerical integration. Concerning the impact of the chain length N we refer to Ref. [13]. We 
determine the maximal can be obtained within the time interval T ~ 

[0, . . . , 4A^/ max( Jo, Ji)]. The results depicted in Fig. (TJa) reveal that at — 2ho the maximal 
concurrence assumes a value close to unity during this time interval and is significantly larger 
than for other frequencies. We refer to this resonance condition as entanglement resonance, 
which holds independently of the other parameters. On the contrary, we find that the height 
and width of the resonance peak depend on the intensities Jq and Ji and on the chain length 
N (not shown in the figure, see Refs. [13]). We also notice the existence of a secondary, smaller 
peak at the driving frequency ojd = ho. In contrast to the main peak, its amplitude strongly 
decreases with the coupling intensity and increasing chain length [T?| . 
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In order to understand this phenomenon quantitatively, we replace the Hamiltonian ^ 
by its time average within a rotating- wave approximation (RWA). At the resonance condition 
ujd = 2ho, we obtain 

J N-l 

= y XI [^n '^n+l + + ^-C , (7) 



with the effective anisotropy 7 = 7J1/2J0. This means that for resonant driving, the time- 
dependent XF-model ^ can be mapped to the static XF-modcl (O without Zeeman fields. In 
both cases, the entanglement generated between the two ending spins is maximal and controlled 
by the parameter 7 and the chain length N. Note that Jo — > corresponds to the infinitely 
anisotropic limit 7 — >■ 00. One important point is worth being mentioned: the concurrence 
approaches unity in the limit of vanishing Jq, i.e. for infinite 7. In this limit, the amount of 
entanglement no longer depends on Ji and 7. An argumentation given in Ref. |13j explains why 
the limit Jg is optimal for entanglement generation. As a rule of thumb, the concurrence 
becomes maximal when minimizing the swap terms in the Hamiltonian ([3]). Let us finally 
mention that in Ref. [T^ , the entanglement scalability with the chain length and the arrival 
time for this entanglement was studied. 



3.2 Modulation of the local fields 

In the following, we investigate the response to driving with local fields. We assume that at 
t = the interaction is switched to a constant value with Jo > 0, Ji = 0. However, the local 
Zeeman field has now both ac and dc contributions, h(t) = ho + hi sin(a;t). Again we assume a 
homogeneous chain, ej = 1 for all j. Proceeding exactly as above, we compute the maximum 
concurrence reachable within the time interval [0, . . . ,4A'^/ max(Jo, Ji)]. Our numerical results 
are drawn in Fig. [ijb) . Also in this case we find an entanglement resonance for = h^ and for 
uJd = 2/10, where the peak at uj^ ~ 2ho dominates clearly. In order to understand this behavior 
we recall the Hamiltonian in the interaction picture, Eq. ([3]), simplified by setting Ji = in 
this case, 

N-l 

^ = Y E + + 7(c'^(*V+a++, + e-'^Wa-a-^,)] (8) 

with S{t) = 2 ds{ho + hi sm{ujdt)). Notice that the driving no longer affects the swapping 
terms. On the other hand, the entanglement in terms of the concurrence is favored by the terms 



+ 

71+1 



-h.c, as pointed out after Eq. ([3]). For this reason, the resonances should assume values 
for which the time-dependent coefficients exp[iZ'(i)] are not averaged to zero. To get further 
insight, we perform a Taylor expansion in powers of {hi/uj), yielding 



where = ri!/fc!(n— fc)! are binomial coefficients. With this at hand we obtain the resonance 
condition 

2/io + {n - 2k)ojd = (10) 

To lowest order. Eq. ^TU\\ is fulfilled for n = 1, which yields oJd = 2/io- This again confirms our 
numerical results depicted in Fig.[T] The next order at 71 = 2 yields the second resonance condi- 
tion LOd ~ ho, which explains the second peak in the figure. Peaks corresponding to resonances 
of higher orders (n > 2) are suppressed for hijujd ^ 1. 

Finally, we emphasize that the swapping terms cr^cr^+i are unaffected by driving the local 
fields. Thus, less concurrence is obtained as compared to the case with time-dependent interac- 
tion described in Sec. 13.11 This feature becomes apparent when comparing the Hamiltonians ([7]) 
and (O or, equivalently. Figs. [TJa) and[TJb). 
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Fig. 1. Maximal concurrence CJjv'^ between the end spins as a function of the driving frequency 
for various anisotropies. (a) Entanglement creation by modulation of the interaction J{t) with the 
parameters hi — 0, Jo — 0, and Ji — O.l/io, cf. Eq. ([2]). (b) Corresponding effect due to modulation of 
the Zeeman fields hn{t), for hi — O.lho, Jo = O.l/io, and Ji = 0. The chain length is A'^ = 6 in both 
cases. 



4 Directing of entanglement current via ac-control 

In Ref. [1] the intriguing possibility of entanglement steering in an isotropic XY chain (7 ~ 0) 
was explored. It was shown how to make excitations in the spin chain propagate into a specific 
direction with a ratchet-like profile of the local Zeeman fields [231 plus a superimposed oscil- 
lation. This suggests the fabrication of a "quantum router" . In this way, distant parties could 
be brought to share a high amount of entanglement, provided that dccohcrcncc is sufficiently 
small. As a further option, the excitation may be split into two parts, each propagating to one 
end of the chain. This results in entangling three distant parties rather than two. In Fig. [2] 
wc sketch a model for the simplest router possible, where the entanglement has two possibili- 
ties of propagation. Corresponding numerical results about splitting of entanglement into two 
directions are depicted in Fig. [31 

To extend and generalize this elementary sketch of a quantum router, we point out how 
multipartite entanglement could be created between many distant parties. For this purpose, we 
suggest choosing a setup in which the entanglement signal has been sent through a driven spin 
chain, split in two. Both tail ends of the chain get thus entangled with the sender of the packet. 
The initial state |100) having been sent by Alice, the final state becomes 

[Alice, Bob, Charlie) = ^[|100) + ^(|010) + |001))j , (11) 

and is now shared by Alice, Bob and Charlie (the entries of the state vectors are given in this 
order). State pT|) is a tripartite entangled state with similarity to a H^-state pHI^ . except 
for the weighting factors. It represents the different quantum coherent histories, i.e. outcome 
possibilities of excitation transfer. The component |100) refers to the case in which Alice has 
not sent any excitation through the driven chain. The latter two terms |010) and |001) result 
from splitting an excitation, which has been sent by Alice before. 

If other driven chains are appended to Bob and Charlie, initial excitations can be split 
further. This results in a final penta-partite entangled state, i.e. a state held between five parties, 
which exhibits yet two more terms. However, the share of the distant parties gets diluted by a 
factor I/V2 every time the signal is split, while Alice's fraction never gets diminished. To get 
over this drawback, Alice can send her excitation to another spin array in a different direction. 
As illustrated in Fig. |H this would result in a quadripartite state, where Bob and Charlie are 
entangled with Beatrix and Carol, 

[Beatrix, Carol, Bob, Charlie) = i [[1000) + [0100) + [0010) + [0001)] (12) 

This state is well balanced, and further splittings keep this balance unaffected if performed 
symmetrically on both ends. We notice too that the cardinality of the states is irow even. 
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Charlie 

Fig. 2. Model for the simplest quantum router possible. A green spin chain allows perfect transfer and 
is not driven, while the propagation in the orange chain is controlled by ac fields. 




t[l/J] 

Fig. 3. Entanglement transport in a spin chain structure corresponding to Fig. [2] Concurrence Co,„ 
[Eq. ([S])] between the spin at Alice's site in the horizontal chain (labeled here n = 0) and the spins in 
the vertical chain. The entanglement is split up, arriving both Bob (n = 12) and Charlie (n — 1) after 
a given time. Here, the chain is isotropic (7 = 0), while the Zeeman fields have a ratchet-like globally 
driven shape as described in Ref. [T]. 

whereas it is odd in the above scenario. To exploit advantages of multipartite distant entangled 
states, one could use a setup such as sketched in Fig. [5] 

5 Experimental realizations and decoherence effects 

A highly accurate and stable confinement of ions in traps has been achieved in the last decade, 
together with the possibility of operating in the ground state of the spatial ionic motion P51l?7] . 
Several proposals to use this amazing tools in order to simulate spin chains have been given [H 
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Fig. 4. Setup for the production of quadripartite distant entangled state, where Alice sends here initial 
singlet into two different driven chains which both are in operated in the splitting mode. 
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Fig. 5. Setup for the production of multipartite (sixteen parties) distant entangled states where Alice's 
singlet is split several times. As above the gray boxes mark output terminals. 



[57] . Here, the use of magnetic field gradients and state dependent laser forces, respectively, yield 
an effective spin-spin coupling mediated by the spatial degrees of freedom. The latter is achieved 
through Coulomb repulsion between ions [35] , or a wire- mediated capacitive coupling of the ions 
residing in neighboring traps |29j . These proposals allow for the accurate simulation of many 
spin Hamiltonians, included Eq. ([T]), provided the mean photon numbers (temperatures) with 
respect to the relevant degrees of freedom are low. It must be noted that in [i] spins as such are 
used indeed, whereas in |37j they are simulated by means of two internal electronic ion levels. 
The latter proposal has been experimentally realised [35] . Furthermore, the first two Fock levels 
of the axial motion can be used to implement spins for the purpose of quantum routing [I], 
given that the ions in the chain can be prepared in the ground state with a probability of 
99.9% [ig. 

In the following, we take into account decoherence, which is the main obstacle for the 
creation of entanglement through resonant driving. It originates from weak but unavoidable 
coupling of the spin chain to a typically large number of uncontrollable external degrees of 
freedom. Irrespective of its physical nature, the environment is usually modeled as a bath of 
harmonic oscillators, each coupling to a system coordinate via its position operator [50^55] . 
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Fig. 6. Maximal concurrence CJ^" for resonant driving [Eq. as a function of the driving frequency 
for various decoherence rates. The anisotropy parameter is 7 = 5, while the chain length is N = 6. (a) 
Modulation of the interaction J(t) with the parameters hi = 0, Jo — 0, and Ji = O.l/io- (b) Modulation 
of the fields hn{t), where hi = O.l/io, ^0 ~ 0-lho, and Ji = 0. 

Here we assume that each qubit n undergoes pure and independent dephasing, i. e. that its 
coordinate crfj couples to a separate bath. 

For weak dissipation, one can eliminate the bath within second-order perturbation theory 
and derive a Bloch-Redfield master equation for the reduced density operator p of the qubits j34j . 
Considering only phase noise on the traps, the equation takes a Lindblad form |35] . 

A ^ 

p^-l_[H,p] + -Y, i<P< - P) , (13) 

where we have assumed that the effective decoherence rate A is the same for all qubits. The above 
master equation has been obtained perturbatively in the dissipation strength. However, it turns 
out to be valid beyond the perturbative regime and therefore provides an accurate description 
of the dephasing dynamics |36| . The heating rate, which would include relaxation terms in 
the master equation, is typically several orders of magnitude smaller than the dephasing rate. 
This justifies the modelling of decoherence by pure phase noise. We use a rather conservative 
estimate and assume A ~ 10~^ J. 

Now we consider the impact of decoherence in the resonance phenomena discussed in Sec.[3l 
In Fig.|n]we plot the concurrence versus the frequency at different decoherence rates. Obviously, 
the concurrence becomes smaller with increasing decoherence. This reduction is related to the 
quantum master equation (jl3p with respect to the time it takes for the entanglement signal to 
reach the ends of the chain. Even so, the entanglement is well preserved for a decoherence rate 
A = O.OOIJ. This is quite a realistic estimate for most physical systems that are suitable for 
experimental verification of entanglement resonance. For this reason, it should be possible to 
observe entanglement transport in an experiment as proposed, even though it will be limited 
by decoherence. 

6 Conclusions 

The transfer of methods established in condensed matter physics to the field of quantum in- 
formation, and vice versa, has recently undergone a blossoming and promoting activity j39j . 
With this work, we have reviewed and extended our proposals pill3j for generating and rout- 
ing entanglement in a spin chain by utilizing time-dependent control of spin-spin interactions. 
Entanglement routing is achieved by a global manipulation of the chain via an external driv- 
ing fields. This method possesses a salient advantage as compared to the rather difficult local 
addressing of individual spins, since that task requires a delicate tuning of intrinsic system 
parameters. An optimal amount of entanglement is generated at a particular resonance, which 
builds a bridge to the familiar theory of coherent destruction of tunneling [5^^111 HT] . 
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Thus, driving with external fields paves the way to control entanglement dynamics in real- 
istic spin chains. We further elaborated on the role of decoherence thereby incorporating the 
unavoidable influence of environmental degrees of freedom in a realistic scenario. With regard to 
an experimental realization, we have sketched the implementation with Penning traps. Another 
possibility consists in using tailored Josephson junction arrays. 

A further advantage of our time-dependent driving scheme is the substantial decoherence 
suppression by use of ac-fields [32ll42l - [44] . This time-dependent manipulation allows the opti- 
mal control in attaining maximal entangled of distant objects [13| . In addition, Landau-Zener 
dynamics has been proposed as well for entanglement generation protocols [IHISS] . Finally, the 
implementation of optimal quantum gates was suggested |47l . All in all, this justifies the ef- 
forts in employing the toolboxes from both driven quantum systems and solid state correlation 
physics in administrating quantum information processing. 

The authors thank the Deutsche Forschungsgemeinschaft for financial support via the collaborative 
research center SFB 484. F.G. acknowledges support by COQUSYS. D.Z. acknowledges partial support 
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